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Abstract 

Let W be a weight-homogeneous planar polynomial differential system 
with a center. We find an upper bound of the number of limit cycles which 
bifurcate from the period annulus of W under a generic polynomial perturba- 
tion. We apply this result to a particular family of planar polynomial systems 
having a nilpotent center without meromorphic first integral. 
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1 Statement of the results 

The real planar differential system 

i = f(x,y), y = g(x,y) (1) 

is said to be (a, f3, uj)-weight-homogeneous if there exist weights a, (3 > and 
uj G R such that for all p > 0, x, y G R holds f(p a x,p^y) = p" +a f{x,y) and 
g(p a x, p 13 y) = p ul+l3 g(x,y). Equivalently, the system (OQ) is weight-homogeneous if 
the foliation 

: f(x, y) dy - g(x, y) dx = 
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is invariant under the dilatation 

$ : R 2 -> R 2 : (x, y) i-> y) = (p a ac, p' 3 y). (2) 

In w/iai follows we shall suppose that the weight-homogeneous system ([T]) has 
a center at the origin. Then the center is global and the open period annulus 
O = R 2 \ {(0, 0)} is a union of periodic orbits. The center needs to be global due 
to its invariance under the aforementioned dilatation. The purpose of the present 
paper is to study small perturbations of such global centers. 

Consider a one-parameter analytic perturbation 

x = f(x,y) + eQ(x,y,e), y = g(x,y) - e P(x,y,e) (3) 

of ([1]) and suppose that P(x,y,e),Q(x,y,e) are polynomials in x, y of degree d, 
depending analytically on the parameter e. For every compact set K contained in 
the real positive half-axis R^ = {(x, 0) : x > 0}, and for |e| small enough, there 
exists an open interval A D K on which the first return map 

IT £ : A -> R + : x i-> U e {x) (4) 

is well defined and analytic 

U £ (x) = x + e k M k (x) + o(e k ). 

Its fixed points correspond to limit cycles of (jSJ). We note that the function M k (x) 
is defined on the whole half-axis IR+ = {(x, 0) : x > 0} and its zeros counted with 
multiplicity provide an upper bound for the number of fixed points of n e on K. 
It is known that the so called higher order Poincare-Pontryagin-Melnikov function 
Mfc(x) allows an integral representation in terms of iterated path integrals [2J [3], S] 
along the periodic orbits 7(2;) of the system ([T]). If the perturbation is generic, then 
k = 1 (Mi ^ 0). The first result of the paper is the following 

Theorem 1 // the first Poincare-Pontryagin-Melnikov function Mi is not identi- 
cally zero, then the perturbed system (0) has at most (d+ l)(d + 4)/2 — 1 limit cycles 
which tend to periodic orbits as e tends to zero. 

In other words, the cyclicity of the open period annulus O with respect to generic 
(such that Mi ^ 0) perturbations of degree d is at most (d+ l)(d + 4)/2 — 1. This 
bound is certainly not exact, as one can easily check in the case of a linear center. 
The above theorem does not make any claim about the number of limit cycles which 
tend to the origin or to the "infinity". 

In the case when the system ([1]) has a polynomial first integral and is moreover 
Hamiltonian, the computation is straightforward (as we have a well known integral 
formula for Mi), see [HI El EI]- However, in general, a weight-homogeneous system 
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with a global center is neither Hamiltonian, nor it has an analytic or even meromor- 
phic first integral. Under the restrictions that the polynomials /, g have no common 
divisor in M.[x, y] and that the origin has no characteristic directions, a result close 
to our Theorem [1] was recently announced in Theorems A, B, C]. We note that 
our result is different, more general and with a much shorter proof. 

The next question addressed in the paper is: are there non-U amiltonian weight- 
homogeneous polynomial systems with a center ? 

The answer turns out to be positive. Theorem [2] provides a large class of non- 
trivial weight-homogeneous systems with a global center to which Theorem [1] ap- 
plies. An explicit example of such a system is given in (jSj) below, and the exact 
upper bound for the number of the limit cycles under a generic perturbation is given 
in Theorem [3J 

To formulate Theorem[2l let hi, h 2 £ M[x, y) be two (a, (3, 5)-weight-homogeneous 
polynomials of the same weighted degree, that is, hi{p a x,p^y) = p s hi(x,y), 
i — 1, 2, for any x, y £ R and any p £ R + . We recall that a, (3, 5 are real numbers 
with a, (3 > 0. We also assume that hi and hi are such that 

• h 2 {x,y) > 0, V(x,y) £ M 2 

• hi(x,y) = h 2 (x,y) = if and only if (x,y) = (0,0) . 
Let a £ C and 3?cr 7^ and put 

H={hi + ihzY^ - ih 2 f, V = -{ln + ^ 2 ) 1 - <T (/ii - ih 2 ) x -°. 



Theorem 2 The system 

x = H y V, y = -H X V (5) 

is a real polynomial (a, (3, 25 — a — (3) -weight-homogeneous planar differential system 
which has a global center. 

The proof of this Theorem is given in Section [2l Clearly some hypothesis can be 
relaxed. For instance, hi, h 2 need not be polynomials. 

Example Put 



hi(x,y) = x 2n + y, h 2 (x,y) = + 4c| x 2n /2, a = 1 - 1/Vl + 4c, 

where n is a natural number n > 1 and c is a real number with c < —1/4. System 
(J5|) takes the form 

x = y + x 2n , y = 2ncx in - 1 , (6) 
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and it has a global center with a first integral 

H=(y + /i+x 2 ") 1 - 1 ^ 3 ^ + /i_£ 2 ") 1+1 /^\ (7) 
where /i± = (1 ± \fl +4c)/2. 

We remark that system (jfJJ) is a (1, 2n, 2n — l)-weight-homogeneous differential 
system. 

To apply Theorem [T] to (jH]), we consider the following perturbed system 

x = y + x 2n + eQ(x,y,e), y = 2ncx 4n ~ 1 - eP(x,y,e), (8) 

where |e| > is a small parameter, P(x, y, e) and Q(x, y, e) are polynomials in (x, y) 
and depend analytically on e and P(x, y, 0) and Q(x, y, 0) are polynomials of degree 
at most 4n — 1. Let { r fh}h be the continuous family of periodic orbits surrounding 
the center. The first Poincare-Pontryagin-Melnikov function Mi can be written as 
follows 

P(x,y,0)dx + Q(x,y,0)dy 

V(x,y) ' 1 } 

(this formula holds true for the system ([3]) too, where V is an appropriate integrating 
factor). We deduce the following 

Theorem 3 Suppose that M\(h) is not identically zero. Then 

(a) the perturbed system (0) has at most n(3n + 1) — 1 limit cycles which tend to 
period orbits as e tends to zero. 

(b) For suitable polynomials P(x, y, 0), Q(x, y, 0), and for all sufficiently small \e\, 
the perturbed system §8§ has at least n{2n + 1) — 1 limit cycles which tend to 
period orbits. 

The bounds are written in terms of the number n associated to system ([6]) which is 
not its degree. The degree d = An — 1, and hence the upper bound for the number of 
limit cycles is (d + l)(3d + 7)/16 — 1, which is better than the one given in Theorem 
[U due to the symmetries of ([6]) . The lower bound can also be written in terms of the 
degree and it is (d+ l)(d + 3)/8 — 1. It is obtained by a direct study of the function 
Mi. The symmetries of system §5§ imply that the lower bound for the number of 
limit cycles ((d + l)(d + 3)/8 — 1) is strictly lower than the corresponding lower 
bound (d(d + l)/2) for a degree d perturbation of a generic Hamiltonian system of 
degree d [5]. 
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2 Proofs. 

Proof of Theorem^] We shall estimate the number of isolated zeros (counted with 
multiplicity) of the first Poincare- Pontryagin-Melnikov function 1(h) = Mi(h), see 
(ED- Consider first the special case Q(x, y, 0) = and P(x, y, 0) = x l yi with i, j > 
and i + j < d. Denote by T e the real foliation defined by 

T e : —g(x, y) dx + f(x, y) dy + ex % yi dx = (10) 

and let IT e be the corresponding first return map. The dilatation $ defined in (J2]) 
transforms T e to the foliation §*T e 

$*JF £ : — g(x, y) dx + f(x, y) dy + ep e x l y j dx = 0, £ = + (3 — ioc — j(3. (11) 

The first return map Tl £p i of the foliation T ep i = $*J r e is therefore conjugated to 

n e = f 1 on 8 ^o^ (12) 

where 0(x) = 0) is the restriction of $ to the cross-section = {(x, 0) : a; > 
<>}• 

Given a foliation with an associated Poincare map II and a diffeomorphism 
then the Poincare map associated to $* T is conjugated to IT. This is due to the fact 
that the Poincare map of a foliation is defined by its flow, see for instance Theorem 
1 in page 207 of [9], and the flows of the foliations T and $*JF are conjugated, see 
for instance Lemma 11 in page 217 of [10] . 

We have 

U s (x) = x + eI(x) + o(e) (13) 
U epi (x) = x + ep £ I(x) +o(e) (14) 
(f)- 1 oH £pl o(j)(x) = p- a (p a x + ep e I(p a x) + o(e)) (15) 
= x + ep £ - a I(p a x) +o(e). (16) 

Therefore, equating the first order terms in e in (fl3|) and (TX6]) we get 

I(x) = p e - a I(p a x), 

for any positive real numbers p, x. This implies, choosing p = a;™ 1 /", that 

I{x) = 1(1) x {a ~ e)/a = x l ~ {u)+p)/a 1(1) x ( ai+ ? j V a (17) 

In a similar way, if we suppose that Q(x, y, 0) = x l yi and P(x,y,0) = with 
*5 j > and z + j < d, we deduce 

I( x ) =I i-M/«j( 1 ) a .(«(i-W+i))/<. ( 18 ) 
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(this computation is omitted). Finally, taking into consideration the additivity of 
the Poincare-Pontryagin-Melnikov function Q with respect to the monomials of 
P(x, y, 0), Q(x, y, 0) we conclude 

_ x l-{u!+/3)/a c _ x (ai+f3j)/a 

where J = G Z 2 : -1 < i < d, < j < d + 1, < i + j < d} and cy el 

The number (d + l)(<i + 4)/2 is the cardinal of the set J . Obviously the function 
I(x) has at most (d + l)(d + 4)/2 — 1 zeros (counted with multiplicity), on the 
interval (0, oo). This completes the proof of Theorem dJ | 

Remark 4 Clearly, the above bound is exact if and only if a and (3 are not com- 
mensurable. The weights a and (3 of any weight-homogeneous polynomial system 
([!]) with a center at the origin are necessarily commensurable. To show this fact just 
take any monomial of f(x,y) with a nonzero coefficient, a.ijX l yi , and any monomial 
of g(x,y) with a nonzero coefficient: b V jix % yi . Since f{p a x,p^y) = p^ a f(x,y) 
and g(p a x, p 13 y) = p UJ+f3 g(x,y), for any x,y,p G M, we deduce the identities: 
ai + (3j = to + a and ai' + f3f = uj + (3. We subtract them to deduce that 
a(i — z' — l) + (3{j — j 1 + 1) = 0, which gives that a and (3 are commensurable. In 
the case when both f(x,y) and g(x,y) have only one monomial with nonzero coeffi- 
cient (f(x, y) = a^yi , g(x, y) = bi'j'X 1 yi ) and such that i — %' + 1 and j' = j + 1 
the weights a and [3 are not commensurable, but the origin is a linear node instead 
of a center. 

REMARK 5 The same result directly follows from the formula (Q for 1(h); it suffices 
to note that H, H x , H y , V are weight-homogeneous functions of appropriate degree. 

Proof of Theorem® The weight-homogeneous degree of the system (a, (3,25 — 
a — (3) follows from straightforward computations. 

The condition that h2(x,y) > implies that the variation of the argument of 
hi(x,y) + ih 2 (x,y) along any closed path I C M. 2 \ {(0,0)} is zero. Therefore for 
every fixed o~ G C the function (hi + ih 2 ) cr has a single valued analytic continua- 
tion on M 2 \ {(0,0)}. From now on we fix some determination of (hx + ih 2 ) a . We 
note that functions H, V defined above, as well the associated differential system, 
do not depend on this particular determination. We may suppose without loss of 
generality that > (otherwise we just replace H by 1/H). Then H has a 
continuous limit at (0, 0) and we may put H(0, 0) = 0. We claim that each level 
set {(x, y) : H(x,y) = e}, e > 0, is a smooth closed curve containing the origin. 
Indeed, the restriction of if on a half-line I starting at the origin is again a positive 
weight-homogeneous function. It follows that {(x, y) : H(x,y) = e} C\ I consists 
of a single point and therefore {(x, y) : H(x,y) = e} is a closed curve. Suppose 
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that dH(xo, y$) = 0. Then the differential of H is zero at any point belonging to 
the half line Iq starting at the origin and containing (xo>2/o)- It follows that H\l Q is 
a constant and moreover this constant equals to = H(0, 0) which is impossible. 
Therefore the level set {(x, y) : H(x, y) = e}, We > 0, is a closed periodic orbit. The 
system has a global center. | 

Proof of Theorem^ According to Theorem[2]the origin of system (jH]) is a center. 
As noted in Remark [5J the first integral H and the inverse integrating factor V are 
weight-homogeneous. More precisely 

Lemma 6 Given any p G C — {0} we have: 

V (px, p 2n y) = p 4n V(x,y), H (px, p 2n y) = p 4n H(x,y), 
for any (x,y) G R 2 - {(0,0)}. 

The proof is a straightforward computation. 
We define the following functions: 

f x l y^ f x l y^ 

Iij(h) = (f> — -dx, Jij(h) = (f) T77 - — dy, 



H =h V(x,y) ' ' J H=h V(x,y) 

where i,j are non-negative integer numbers. By the expression of 1(h) and taking 
into account that P(x, y) and Q(x, y) are polynomials with real coefficients of degree 
at most 4n — 1, we have that: 

(h) + PijJijih), (19) 

0<i+j<4n-l 

with a,ij and real numbers and i,j are nonnegative integer numbers. 
Lemma 7 Given any p, h G K with h > and p ^ 0, we have: 

i id ( P in h) = | P r +i+2 ^- 4 %-(/i), j tj ( P in h) = ipi^^-^j^h), 

where \p\ stands for the absolute value of p. 

Proof. The change of variables x i— ► \p\x, y i— > p 2n y in the integrals Jy and J^- 
(which preserves the orientation of the oval {if = h}) implies 

J..(n^ h ) = I dx = I I |i+l+2nj-4n , 

MP j lH^ h V(x,y) dX f H J Pl V(x,y) dX 
= \p\ l+1+2n] - An l,j(h). 

X ' yJ A„. - I I \i+2n(j+l)-4n X% V 3 



JH= P ^ h V(x,y) J H=h V(x,y) 

= |pr +2n(j+1) - 4n J 4J (/i). 
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From the previous lemma we give the form of the functions hj(h) and Jij(h). 
The same argument given in the proof of Theorem [1] holds: we choose h = 1 and 
p = in the expressions given in Lemma [3 

Lemma 8 Given any h > we have: 

i+l + 2nj-in i + 2n(j + l)-4,n 

Jy(ft) = k 4n Jij{h) = h 4n Jy(l). 



We get that hj(h) and Jij(h) are monomials of h up to a fractional power. Hence, 
we only need to determine how many of these functions are linearly independent, 
so as to know how many zeroes can have 1(h). We first determine which of the 
functions hj(h) and Jij(h) are identically zero. 

Lemma 9 If i is odd, then hj(h) =0. If i is even, then Jij(h) = 0. 

Proof. Let us consider the change of variables x \— > — x and y i— > y in the integrals Iij 
and Jij. We note that this change of coordinates reverses the orientation of the oval 
{H = h}. We denote the oval with reversed orientation by — {H = h}. Moreover, 
from Lemma [6] we have that this change of coordinates leaves V(x, y) and H (x, y) 
invariant. 



Taking i odd, we deduce that hj(h) needs to be identically zero, and the same is 
true for Jij(h) taking i even. | 
Lemma [9] can also be proved using Green's Theorem and analogous reasonings. 

We are going to characterize some of the functions hj(h) and Jij(h) which are 
not zero at any point of h > 0. 

Lemma 10 For any k, £ non-negative integers we have: 

(h) j£ 0, J 2k+ iA h ) £ o, 

where we recall that 

I W (h) : T{ ^y dx, J 2M Mh) := T{ - y dy, 

with H(x, y) as given in (J7J) and V(x, y) := y 2 + x 2n y — cx 4n . 
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Proof. We note that the orientation of system (jS]) over the oval H = h is clockwise 
and that V(x, y) is strictly positive over all the oval H = h. 

Let us denote by x\{h) and x 2 (h) the intersections of the oval H = h with the 
horizontal axis (y = 0) with x\(h) < and x 2 (h) > 0. We denote by {H = h} y<0 
the half part of the oval below the horizontal axis oriented from X2(h) to x\{h) and 
by {H = h} y>0 the half part of the oval above the horizontal axis oriented from 
Xi(h) to x 2 (h). We have: 

r x 2ky2l+l rxi(h) x 2ky2£+l rx 2 {h) % 2k 21+1 

hk,2£+i{h) = f —, r-dx = / — — -dx + / — — -dx 



H =h V{x,y) J X2(h) V{x,y) J Xl{h) V{x,y) 

xa{.h) x 2k y2£+l rx 2 (h) % 2k 21+1 

dx + — ; — dx 



i Xl (h) V{x,y) ' ' J Xl(h) V(x,y) 

where the first integral is done over the path {H = h} y<0 and the second integral 
is done over the path {H = h} y>0 . Therefore, the first integral is strictly negative 
and the second integral is strictly positive and we are adding two positive values, 
due to the minus sign. Hence, l2k,u+i(h) > and it cannot be zero at any point. 

Analogously for J2k+i,2e(h), we define yi(h) and yzih) the intersections of the 
oval H = h with the vertical axis (x = 0) with yi(h) < and y 2 {h) > 0. We denote 
by {H = h} x< Q the half part of the oval at the left of the vertical axis oriented 
from yi(h) to yz(h) and by {H = h} x>0 the half part of the oval at the right of the 
vertical axis oriented from y 2 {h) to yx(h), and we have: 

r x 2k+ly2i ryi(h) x 2k+ly2£ ryn(h) x 2k+l y 2l 

J2k+i,2t{h) = <b — ( r dy = — dy + — dy 

J H =h V{x,y) J y2{h) V{x,y) J yi(h) V(x,y) 

V2(h) 2k+l -.21 r y 2 (h) 2k+l -.21 

^ y ^y _|_ / ^ y ^y 

hx{h) y(x,y) J yi (h) V{x,y) 

where the first integral is done over the path {H = h} x> o and the second integral 
is done over the path {H = h} x<0 . Therefore, the first integral is strictly positive 
and the second integral is strictly negative and we are adding two negative values, 
due to the minus sign. Hence J 2 k+i,2i{h) < and it cannot be zero at any point. | 

In fact, we have been able to numerically prove that the integrals I 2 k,2i{h) ^ 
and J 2 k+i,2t+i{h) ^ for some particular fixed values of the integers k, £, with 
k > and £ > 0, in some fixed cases of the function H(x, y) defined in (j7j). To show 
this fact, we have parameterized the oval H(x,y) = h by (x + (r),y(r)) when x > 
and (se_(t), v(t)) when x < 0, with 

x±(t) = ±h*^ (t + (r + fi-)~^, y(r) = hh r (r + (r + 

where fi± = (1 ± y/1 + 4c)/2, <j — 1 — + 4c and the rank of the parameter r 
is all the real line r £ (— oo, +oo) in both parts of the oval. When we write or 
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hh we mean the positive real root. 



Next lemma shows that the possible nonzero values of the integrals Jij(h) are 
redundant, since there is an integral of the form Iij(h) which corresponds to the 
same monomial. 

Lemma 11 The first Melnikov function can be expressed by: 



where G M. and I = {(«, j) 6 Z 2 : < i, j < An — 1, % is even, i + j < An — 1}. 

Proof. The expression of 1(h) given in ffT9l) ensures that this function is a linear com- 
bination of the integrals hj(h) and Jij{h). We are going to show that any possible 
monomial expressed by a Jij{h) can also be got by a monomial of Ii>j>(h). Lemma 
M gives that only the integrals Jij(h) with i odd need to be considered, hence we 
take any two nonnegative integers such that i is odd and % + j < 4n — 1. We 
define %' = % — 1 and j' = j + 1 and we have that both i' and j' are nonnegative 
integers strictly lower than An and i' + j' — i + j < An — 1. Moreover, is even in 
accordance with Lemma [9] and (i',f) G X. Hence, any monomial given by a Jij(h) 
is also expressed by a /^-(/i) with z even. | 

To end with the proof of Theorem [31 we first count the cardinal of X and we are 
going to show that jjX = 2n(2n + 1). We take any even nonnegative integer i from 
to An — 2, that is, we take 2n possible values of k with i = 2k and given a fixed 
< k < 2n — 1 we can take any value of j from to An — 1 — 2k. Hence, 

2n-l 4n-l-2fc 2n-l 



We note that the cardinal of X is an upper bound for the number of independent 
monomials given by the nonzero Iij(h) because it may happen that two elements 
of Iij{h) give rise to the same monomial, that is, it may happen that (i, j) G X, 
{i',j') G X with (i, j) ^ but i + 2nj — il + 2nj'. To end with the proof of 

part (a) of Theorem [3l we need to characterize the number of repeated exponents 
corresponding to different indexes (i,j) G X, (i',j') G X. We note that if two such 
pairs give i + 2nj = i! + 2nj', then j and j' can differ at most by one, because if 
they differ by two or more then i and i' differ by An or more which is not possible 
since < i, i' < An — 1. We assume that j' = j + 1 and we have that %' — i — 2n. 
Hence, fixing (i, j), the condition to have a repeated exponent is that i > 2n and 
we already have that j < An — 1 because i + j < An — 1 . Let us count how many 



(m')gx 




1 = (4n - 2fc) = 2 + 4 + 6 + . . . + An = 2n{2n + 1). 



k=0 
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of these indexes we do have in X. We fix k such that i = 2k and k goes from n to 
2n — 1. Given such a A; we can take, as before, any value of j from to 4n — 1 — 2k. 
Therefore, the number of repeated exponents is: 

2n-14n-l-2fc 2n-l 

1 = ^(4n - 2k) = 2 + 4 + . . . + 2n = n(n + 1). 

k=n j=0 k—n 

We conclude that the number of different exponents associated to the indexes of X 
is 2n(2n + 1) - n(n + 1) = n(3n + 1). 

We have given an upper bound for the number of independent functions in which 
1(h) can be split as a linear combination. Since this upper bound is n(3n + 1), we 
have that an upper bound for the number of isolated zeroes of 1(h) is n(3n + 1) — 1. 

The proof of part (b) in Theorem [3] comes from Lemma [TTJ], which ensures that 
only half of the 2n(2n + 1) functions Iij(h) are ensured to be different from zero, 
the ones with an even % and an odd j. To end with, we only need to see that two 
such functions give place to two independent monomials. Let (2k, 2£ + 1) G X and 
(2k', 21' + 1) G X with (2k, 2£ + 1) ^ (2k', 2£' + 1) be such that 2k + 2n(2£ + 1) = 
2k' + 2n(2£' + 1) and we will get a contradiction. If k = k' then I = £', so we can 
assume that k < k'. We have 2n(£ — £') = k' — k which gives that £ > £' + 1. This 
inequality gives that k' > 2n + k and since k > 0, we have that k' > 2n, which is 
impossible since 2k' < An — 1. Hence, any two exponents given by different functions 
hk,2e+i(h) are independent. Choosing adequate parameters in system (Q we get that 
1(h) can always be given as a sum of at least n(2n + 1) independent functions and, 
adapting parameters, it can always be chosen with at least n(2n + 1) — 1 different 
isolated zeroes. | 



A The Andreev theorem 

We already noted that Theorem [2] implies that the origin of the system ([6]) is a global 
center. It is interesting to note that the same result follows also from a classical 
theorem of Andreev[TJ which we briefly illustrate. This approach can be useful in 
other cases of nilpotent centers, which are not covered by Theorem [2j 

Theorem 12 (Andreev's Theorem [lj) Let F(x,y) andG(x,y) be analytic func- 
tions in a neighborhood of the origin of order > 2 and let the origin be an isolated 
singularity of the differential system: x = y + F(x,y), y = G(x,y). Let y = <p(x) 
be the solution of the equation y + F(x,y) = such that (f)(0) = 0. We denote by 
= G(x, 4>(x)) and A(x) = div(x, 4>(x)) and we develop them in a neighborhood 
of x = 0: 

£(x) = a lX kl + O (x kl+1 ) , A(x) = a 2 x k2 + O (x k2+1 ) , 
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where a% 7^ 0, k\ > 2 and a 2 7^ 0, k 2 > 1 or A(x) = 0. TTie origin of the differential 
system is monodromic if, and only if, ct\ < 0, k\ is an odd number, and one of the 
following three conditions holds: 

(a) k\ = 2k 2 + l and c? 2 + Aa^^ + 1) < 0, 

(b) ki < 2k 2 + 1, 

(c) A(x) = 0. 

The function associated to (jSJ) is £(x) = 2ncx An ~ l , so cti = 2nc and k\ = 
An — 1. The divergence of system ([S]) is div(a;, y) = 2nx 2n ~ 1 , so A(x) = 2nx 2n_1 and 
we have «2 = 2n and k 2 = 2n — 1. Since n is a natural number and c < — 1/4, we 
already have that a% < and k\ is an odd integer. The condition (a) of Theorem 
[T2]is satisfied because 2k 2 + l = An-1 = h and al + Aai{k 2 + 1) = 4n 2 (l + 4c) < 0. 
We conclude that the origin of system is monodromic. Finally, the origin of 
the system ([6l) is time-reversible because it is invariant by the change (x, y, t) — ► 
(—x,y, —t), which implies that it is a center. Moreover the center is global, as the 
system is weight-homogeneous. 

Acknowledgment. We thank I.D. Iliev for the critical remarks. 
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